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Anisotropic Priedmann-Robertson- Walker universe from nonlinear massive gravity 
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Motivated by the recent no-go result of homogeneous and isotropic solutions in the nonlinear 
massive gravity, we study fixed points of evolution equations for a Bianchi type-I universe. We find 
a new attractor solution with non-vanishing anisotropy, on which the physical metric is isotropic 
but the Stiickelberg configuration is anisotropic. As a result, at the background level, the solution 
describes a homogeneous and isotropic universe, while a statistical anisotropy is expected from 
perturbations, suppressed by smallness of the graviton mass. 



Introduction. General relativity, which describes 
long range gravitational interactions, is in agreement 
with current experimental and observational data. On 
the other hand, from a theorist's point of view, it is an 
interesting question whether the range of gravity can be 
consistently made to be finite, or equivalently, whether 
graviton can have a non-zero mass. This question has 
also phenomenological relevance, since the graviton mass 
may introduce new terms that mimic dark energy and 
thus may source the late-time acceleration of the uni- 
verse. 

Assigning a mass to graviton has been one of the most 
challenging problems in classical field theory for the last 
70 years. The linear theory of Fierz and Pauli [l| gives 
rise to a discontinuity in the observables Q , which can be 
alleviated by nonlinear terms However, such terms 
generically introduce the so-called Boulwarc-Dcscr (BD) 
ghost [H , spoiling the stability of the theory. 

Only recently, a nonlinear extension of the massive 
gravity theory has been introduced [H, 0] , where the BD 
ghost is systematically removed by construction. The 
theoretical and phenomenological possibilities brought 
this theory a significant attention. 

The nonlinear massive gravity allows self-accelerating 
open Fricdmann-Robertson- Walker (FRW) solutions 
with the Minkowski fiducial metric Q as well as 
flat/closed/open FRW solutions with general FRW fidu- 
cial metric |9|. Unlike the other branch of solutions 
[lol Hl| , these backgrounds evade the Higuchi bound [l2[ 
and thus are free from ghost at the linearized level even 
when the expansion rate is significantly higher than the 
graviton mass. This is because there are only two prop- 
agating modes on these backgrounds. However, these 
constructions exhibit aghost instability at nonlinear or- 
der in perturbations [13l | . This is a consequence of the 
FRW symmetries; in order to obtain a stable solution, 
some of these symmetries need to be broken. 

An inhomogeneous background solution was obtained 
in 0, where the observable universe is approximately 
FRW for a horizon size smaller than the Compton length 
of graviton. Similar solutions with inhomogeneities in 
the Stiickelberg sector, meaning that the physical metric 



and the fiducial metric do not have common isometries 
acting transitively, were found in [ljj]. Note that those 
inhomogeneous solutions cannot be isotropic everywhere 
since isotropy at every point implies homogeneity [lij ]. 
Note also that cosmological perturbations can in prin- 
ciple probe inhomogeneities in the Stuckclbcrg sector. 
For example, generic spherically-symmetric solutions are 
isotropic only when they are observed from the center of 
the universe. 

The goal of the present paper is to introduce an al- 
ternative option, where the assumption of isotropy is 
dropped but homogeneity, i.e. the cosmological princi- 
ple, is kept. In a region with relatively large anisotropy, 
we find an attractor solution. On the attractor, the phys- 
ical metric is still isotropic, and the background geome- 
try is of FRW type. Hence, the thermal history of the 
standard cosmology can be accommodated in this class of 
solutions. However, the Stiickelberg field configuration is 
anisotropic, which may lead to effects at the level of the 
perturbations, suppressed by smallness of the graviton 
mass. 

The action and background. We consider a simple 
description of the universe at present time. We assume 
that the late-time acceleration is sourced by a cosmologi- 
cal constant A, as well as the contribution from the gravi- 
ton mass. (Setting A — corresponds to self-accelerating 
solutions as in the example shown in Fig[T]) For this pur- 
pose, the vacuum configuration is sufficient. The action 
with the graviton mass term, constructed by imposing 
the absence of BD ghost in the decoupling limit [6[ , is 



1 = 2 / d 4 xy^[R - 2 A 



where 



'-(C 2 + a 3 C 3 



(1) 



C 2 = [/C] 2 - [K 2 ] , 
£ 3 = i([/C] 3 -3[/C][/C 2 ]+2[/C 3 ]), 
d = ^([/C] 4 - 6[/C] 2 [/C 2 ] + 3[/C 2 ] 2 + 8[/C][/C 3 ] - 6[/C 4 ]) 



(2) 



In the above, the square brackets denote trace operation 
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and 

K£ = W-(Vr T ff v - (3) 

Here, g^ is the physical metric, while the space-time 
tensor is the fiducial metric, whose vacuum expecta- 
tion value gives rise to the breaking of general coordinate 
invariance. 

For the physical metric, we adopt the axisymmetric 
Bianchi type-I metric, which is the simplest anisotropic 
extension of FRW ansatz 

gf}dx^dx u = -N 2 dt 2 + a 2 [e 4,J dx 2 + er 2 ' J S pq dy l dy 1 }, (4) 

where N, a and a are functions of t, Greek indices de- 
note space-time coordinates, while i, j = 2, 3. As for the 
fiducial metric, we assume the flat FRW form as 

ffiv = -^d^d^+^id^d^+s^d^d^), (5) 

where n and a are functions of the temporal Stiickelberg 
field (j>°. The form ([5]) includes a de Sitter fiducial as a 
special case, with Hf = a/ an = constant. 

Varying the Stiickelberg fields around the background 
value 4> a = x a + 7r a , the variation of the mass term up to 



first order is 

/ = + M 2 Pl m 2 g J d i xNa 3 nir°£ <p + 0[(7r a ) 2 ], (6) 

giving the equation of motion 

£ =J { f (H + 2^-H f e- 2,J X) 

+ 2 Jl y) (if - £ - H f e a X) = 0, (7) 

where 

J^ x) = ll -2 l2 e°X + l3 e 2 °X 2 1 

4 V) ^7i-72 {e- 2 ° + e°) X + 73 e" CT X 2 , (8) 

with 7i = 3 + 3 a 3 + a 4 , 72 = 1 + 2 a 3 + a 4 , 73 = a 3 + a 4 , 
H = £ = and X = — . The expansion rate for 
the fiducial metric Hf is related to the invariants of the 
field space metric, and is independent of the choice of the 
background values of <p a . Thus, Eq.([7]) can be interpreted 
as an algebraic equation for a (or equivalently for X), 
instead of a differential equation. 

Varying the action with respect to g^ u , the field equa- 
tions for the physical metric are obtained as 



where 



3(ii 2 -£ 2 ) - A = 
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1 2 [- (3 7l - 3 72 + 73) + 71 (2 e CT + e~ la ) X - 72 {e 2a + 2 ) X + 73 X' j 
i {e~ 2 a -e°)X [ 7l - 72 (e CT + r)X + 73 re a X 2 ] , 



X 
NX 



(9) 



+ H 



(10) 



Additionally, there is also an equation for H , which can 



r 



be recovered by combining Eq.JT]) with Eq.©. 

Fixed Points. We consider a de Sitter fiducial metric 
(Hf = const.) and seek solutions with H = £ = X = 
0. The constancy of X allows us to express H as H = 
Hf X r. In this setup, the independent equations become 



3A - (3 7 i - 3 72 + 73) + 7i(2e ff + e- 2a )X - [^{2e~ a + e 2a ) + 3r 2 u~ 2 ] X 2 + 73A 3 = 0, (11) 
(e CT - 1) [ 7l - l2 {r + e a )X + l3 e°rX 2 ] = 0, (12) 
71 (3r - 2e CT - er 20 ) - 2 72 [(2e a + e- 2a )r - (e 2a + 2 e -°)] X + 73 [(e 2a + 2e~ a )r - 3] X 2 = 0, (13) 

I 

rewritten as 

(7 1 -7 2 Ae CT )(e CT -r)(re 2CT -l)=0. (14) 
Considering (fT4")) as an algebraic equation for e a , there 



where A = ^—^ and u = are dimensionless parame- 
ters. 

For cj = 0, the set of equations is reduced to that for 
isotropic configurations, which was already investigated 
in M,M- Assuming a 7^ and using (Tj"2"j) . Eq. ([T3")) can be 
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arc three solutions: 



71 -1/2 

li X 



(15) 



We now consider each solution separately. 

Case I. e 17 = —^x- Using this solution in Eq. ([T2"j) 
gives X = 71/72, implying a = 0. In other words, this 
solution is isotropic and thus is not of our interest. 

Case II. e a = r. In this case, Eq. p^l) gives 



(r- 1) [71 -2 l2 rX + l3 {rXf 



0. 



(16) 



This equations have two solutions; The first solution 
is r = 1, and leads to isotropy a = which is not 
our interest. The second solution gives rX = (72 ± 
7i73)/73; which reduces Eq. (fTTj) to a nontrivial 
constraint between the parameters of the theory. Since 
this case requires a fine-tuning of a parameter, it is not 
of our interest either. 
Case III, r = e~ 2c 



reduced to 

lie* 7 - l2 {e 2a + e-°)X + l3 X 2 = 0. (17) 

while Eci. (fTTj) becomes 

(3A - 3 7 i + 3 72 - 73) + 7i (e~ 2a + 2e a )X 
-[ 72 (2e- CT + e 2a ) + 3 e~ Aa yT 2 \X 2 + 73A 3 = 0. (18) 

Combining these two equations, we obtain an expression 
linear in X, 

Y _ 3 7i + [7i72 - 73 + 373(72 ~ 7i + A)]M 2 e 3g 

[e° + e^-) [3 72 + (7 2 2 - 7i73)]M 2 e 3ff ' 1 ' 



and an equation which only depends on a 



c Q + Cl e 3a + c 2 e 6a + c 3 e 9a = 0, 



(20) 



With this solution, Eq. (fT2|) is where 



co = 372 (71 + 37! - 37172 - 7273 + 372A) fj 2 - 97 2 , 

ci = (72 - 7173) [-6(371 - 372 + 73 - 3A) + (jl + 372 - 37172 - 7273 + 372A) fj 2 ] /i 2 , 

C2 = [372 + (272 - 27i7 3 )/i 2 ] (jf + 37! - 37172 - 7273 + 372A) n 2 - (7172 - 37173 + 37273 - 73 + 373A) V 4 , 

c-3 = (72 ~ 7i73) (7? + 3 72 - 37172 - 7273 + 372A) // , (21) 



For a given set of parameters (03, a^, A, fi), one can solve 
the cubic equation ([2"U]) for e 3cr and then use fit?!) to cal- 
culate the corresponding value of X. 

Local Stability. We now introduce homogeneous 
perturbations around the fixed point described in the 
third case above. 

H = H f [r X + ehx{t) + 0(e 2 )], 
a = a +ea 1 (t) + 0(e 2 ), 
X = X + eX 1 (t) + O(e 2 ), (22) 

where (Xq, <jq, tq = e~ 2<T °) is the background represent- 
ing the fixed point. Adopting this expansion, we calcu- 
late the equations of motion up to 0(e). At linear order, 
<j\ can be decoupled from the remaining 0(e) quantities, 
and a second-order evolution equation is obtained as 



+ 3X e- 2a -^+M 2 a 1 = 0, (23) 



where 

2 Xfo 2 e-*°° f di (3d 1 -d 2 )(6 + d 1 n 2 ) \ 
2 V 2d 2 -di 2 M 2 / ' 
rfi = (e 3CT0 -l) [ 7 2 - 73e ao A ] , 

d 2 = (e 3 °«-l) [ l2 (3 + 2e 3a ")-5 l3 e ao X ] , (24) 

while the dimensionlcss time coordinate r is defined by 
dr = HfNdt. The fixed point is locally stable if 

M 2 > 0. (25) 

Global Stability. To study the global stability of 
the fixed point, we consider an example with 

A = 0, fi = 20, a 3 = -l/20, a 4 = 1 , (26) 

for which the local stability condition (|25[) is satisfied. 
For this parameter set, there is only one set of real solu- 
tions to the equations of motion (TiTjl - lfTB"]) 

A~4, e ff ~i, r~4. (27) 

In order to determine the phase flow, we first reduce the 
system of equations. Using Eq.®, the first of Eq.@ and 
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FIG. 1: The phase flow for (a, E) for parameters (03, Q4, 
A, /1) =(—0.05, 1, 0, 20). The flow is directed toward the red 
dot at (a, E) = (0.5, 0), which is the fixed point obtained by 
solving Eas.l|TT )l -(|T5 1l . 



their time derivatives, we can express X, H and their 
derivatives in terms of a and S. Since these equations 
are nonlinear, there are always more than one solution. 
For the parameter set (|26|) . we find that there are three 
branches of solutions which give X > and H > 0. 
For each branch, we use this solution with the second of 
Eq.©. As a result we obtain, for a set of (a, £), the cor- 
responding set of (it, S) pair. Out of the three branches, 
only one contains an attractor. The flow corresponding 
to this branch is shown in FiglTJ The flow proceeds to- 
wards the fixed point we found in Eq. (j2"T)) . 

Discussion. The recently introduced nonlinear mas- 
sive gravity theory [5, 6] provides a new framework to ad- 
dress some of the intriguing issues in theoretical physics 
and cosmology, such as the possibility of finite range grav- 
ity and the mystery of dark energy. Although the theory 
admits homogeneous and isotropic solutions, these suffer 
from an unavoidable nonlinear ghost [13[ or a linear ghost 
fioj . Since this is a consequence of the FRW symmetries, 
either homogeneity or isotropy needs to be broken in or- 
der to obtain a stable solution. 



In the present paper, for the first time in the non- 
linear massive gravity theory, we explored regions with 
relatively large anisotropy for homogeneous attractor so- 
lutions. The classification of fixed points revealed the 
existence of a single anisotropic attractor. The local and 
global stability analyses indicate that, a universe with a 
sufficiently large anisotropy at the onset of the late-time 
accelerated expansion should flow to this point. 



A very interesting implication is that the scale fac- 
tors corresponding to the two directions differ only by a 
constant normalization, thus the expansion rate is com- 
pletely isotropic. In general relativity, such a solution 
will be identical to an isotropic universe, up to a co- 
ordinate redefinition. Conversely, in nonlinear massive 
gravity, such a redefinition cannot remove the anisotropy 
completely; it is merely shifted from the physical metric 
to the fiducial metric. Although the background met- 
ric is of FRW type, the signature from the anisotropy 
will be imprinted on the spectrum of cosmological per- 
turbations. The statistical anisotropy signal is expected 
to be suppressed by smallness of the graviton mass m g . 
The type of anisotropy, i.e. statistical anisotropy for per- 
turbations without anisotropic background expansion, is 
totally new. For example, none of the anisotropic infla- 
tion scenarios [l6| has this type of anisotropy. Detailed 
analysis of perturbations and comparison with observa- 
tional data are worthwhile pursuing. As the first step, a 
preliminary analysis of perturbations indicates that the 
anisotropic attractor solutions found in this paper are 
free from ghost for a range of parameters 17 1. 

As a sensible effective theory, the theory can be reliable 
only below a cutoff scale A„ = (Mpi m^ 1 ) 1 / 71 , which is 
much larger than the graviton mass m g for n > 0. The 
exact form of this scale will be determined by a detailed 
analysis of perturbations in a future publication [Tt| ■ For 
the scenario discussed in the present work, we associated 
the late time acceleration with the graviton mass term, 
thus the present expansion rate is taken to be of the 
order of the mass term. For instance, for the case given 
in Eqs. (|26H27[) . the constant expansion rate on the fixed 
point is H = 0.8 m g . This is clearly well below the order 
of the cutoff scale A„ for n > 0. 
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